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Static and dynamic testing of the NASA 0.7 percent scale Hybrid Wing Body (HWB) 
configuration was conducted in the Rolling Hills Research Corporation water tunnel to 
investigate aerodynamic behavior over a large range of angle-of-attack and to develop 
models that can predict aircraft response in nonlinear unsteady flight regimes. This paper 
reports primarily on the longitudinal axis results. Flow visualization tests were also 
performed. These tests provide additional static data and new dynamic data that 
complement tests conducted at NASA Langley 14- by 22-Foot Subsonic Tunnel. HWB was 
developed to support the NASA Environmentally Responsible Aviation Project goals of 
lower noise, emissions, and fuel burn. This study also supports the NASA Aviation Safety 
Program efforts to model and control advanced transport configurations in loss-of-control 
conditions. 
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Fourier coefficients 
deficiency function parameters 
wing span, ft 

normal and side force coefficients 
rolling and yawing-moment coefficients 
pitching -moment coefficients 
mean aerodynamic chord, ft 
deficiency functions 

frequency, Hz 

reduced/, nbf IV , or ncf IV 

no. of harmonics in Fourier expansion 

number of data points 

roll, pitch, and yaw rates, rad/sec 

coefficient of determination 

reference area, ft^ 
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time, sec 
velocity, fps 

angle-of-attack, rad or deg 
initial angle-of-attack in forced 
oscillation experiments, rad or deg 
sideslip angle, rad or deg 
state variable 
standard error 
dummy integration variable 

non-dimensional time constant, — 

h 

angular frequency, rad/sec 
in-phase coefficients, aor P 

out-of-phase coefficients, p, q, or r 

residual and total sum of squares 



Subscripts 
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a 
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amplitude 

aero forces and moments: N, Y, I, m, or n 
measured value 


Superscripts 

^ = estimated value 

= mean value 
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I. Introduction 

T his study contributes to development of NASA Hybrid Wing Body (HWB) aerodynamic models and to the 
investigation and development of mathematical models that can predict aircraft response in nonlinear unsteady 
flight regimes. A nonproprietary version of the HWB (HWB-N2A) is used in this study. Data from the tests provide 
additional information that complement tests conducted at NASA Langley 14x22 wind tunnel facility, for the NASA 
Environmentally Responsible Aviation (ERA) Project \ ERA Project goals are to develop advanced configurations 
with lower noise, emissions, and fuel burn^. This study also supports the NASA Aviation Safety Program, Vehicle 
Systems and Safety Technology Project, where project goals include model and control law development for 
advanced transport configurations in loss-of-control conditions' "^. 

Application of a more general modeling approach to the HWB is presented in this paper. This approach allows 
for the inclusion of nonlinear behaviors and unsteady behaviors. The method takes advantage of indicial models that 
have been successfully applied to both transports and fighter aircraft over many years^'^^. An advantage of the 
method is that the models retain conventional static and dynamic terms but replace conventional acceleration terms 
with indicial functions. This mathematical model structure allows unsteady aerodynamic behaviors to be represented 
under conditions that conventional series expansions of aerodynamic forces and moments cannot model. High 
fidelity aerodynamic models support development and application of simulation, analysis, and control design 
technology to enhance safety of flight of advanced transport configurations. 

Static and dynamic testing of the NASA 0.7 percent scale Hybrid Wing Body (HWB) configuration was 
conducted in the Rolling Hills Research Corporation (RHRC) water tunnel to investigate aerodynamic behavior over 
a large flight regime and to develop models that can predict aircraft response in nonlinear unsteady flight regimes 
A very fine grid test matrix was used to fully characterize the aerodynamics. Elow visualization tests were also 
performed using dye ports along the leading edge of the wings and using a movable wand to investigate flows over 
the body. This paper presents static and dynamic test results, flow visualization graphics, and modeling analysis for 
the HWB. Data and analysis are primarily provided for the longitudinal axis in this paper. An additional study 
including comparisons of wind tunnel and water tunnel results is provided in Ref. 17. 

II. Water Tunnel Model 

The HWB model was designed to perform experiments in the RHRC research water tunnel. Model 2436. In 
order to maintain a wingspan equal to 75% of the tunnel width, a model scale of 0.007 to 1 was used. This resulted 
in a model wingspan of approximately 18 inches. The model is sting mounted, with the balance center located at the 
HWB moment reference center, which is located at the 38.7% mean aerodynamic chord. The balance is also located 
at buttline and waterline 0.0. Eigure 1 shows three views of the general layout of the sting and balance inside the 
model. Reference parameters are: S = 70.263 in^, h = 17.892 in, and c = 7.398 in. 

The model was constructed using several different materials. The main parts are shown in the CAD drawings as 
different colors. The main fuselage was made using a rigid plastic in a stereolithography (SLA) process. This was 
desirable because it allowed relatively easy and inexpensive routing of the dye tubes on the forebody. The nacelles, 
tails and fuselage cover plate were made of SLA material. These parts were painted and sealed with clearcoat to 
prevent water absorption. Because the wings and tails are so thin at this model scale, warpage was a major concern. 
Eor this reason, the wings were initially machined with a computer controlled milling machine out of a stiff nylon 
material and later rebuilt from solid aluminum and then hard anodized. The paint was removed from all of the other 
parts and they were primed with a marine barrier coat and then painted with marine gel coat. 

The balance mounting plate, which was originally a flat surface built into the SLA fuselage, was not rigid 
enough and allowed the model angle-of-attack to change with relatively small pitching moment loads. To fix this, 
the surface was machined out of the model and a new stainless steel balance mounting plate was machined to fit in 
the pocket. Eigure 1 shows the mounting plate in a salmon color at the front of the balance. The mount was epoxied 
into place in the body, which eliminated the deflections under load. The basic water tunnel model performed very 
well for many hours of continuous use. 


III. Measurements 

A series of static tests ^ on the HWB were completed in 2012 at the NASA Langley 14-by 22-Eoot Subsonic 
Tunnel. The water tunnel was used in 2013 to obtain some early dynamic data including some comparison static 
data and flow visualization data. Measurements were made for all three axes with an emphasis on the pitch-axis 
forced oscillation tests. Tests in the water tunnel were conducted at a speed of 11 in/sec, corresponding to a 
Reynolds Number of 5.8x10"^ based on mean aerodynamic chord. The steady data were obtained for cruise and 
droop- wing configurations over a wide range of angle-of-attack from -10° to 30°, and sideslip angles from -10° to 
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10°. The dynamic tests included one degree-of-freedom forced oscillation around the pitch, roll, and yaw axes at 
different nominal values of the angle-of-attack, for different frequencies and amplitudes. Commanded oscillations 
were in the form of single sinusoids, Schroeder sweeps, and ramp-and-hold inputs. Dynamic test results presented in 
this paper are for the cruise configuration under sinusoidal forced oscillation. 

The longitudinal oscillatory data were taken for the cruise configuration, at ten nominal values of angle-of-attack 
in the range 0° to 20°, three amplitudes of 2.5°, 5°, and 10°, and eight non-dimensional frequencies ranging from 
0.05 to 0.30 or dimensional frequencies in the range of 0.0240 Hz to 0.1437 Hz. During the test runs, data were 
obtained for approximately 30 oscillation cycles. Oscillatory data were sampled at 15 Hz with 5 Hz analog and 2Hz 
digital filters. The resulting data were further filtered with a 0.6 Hz low pass digital filter to remove unwanted 
frequency content. The filter was run in both directions to ensure no phase error was added to the data. Initial cycles 
were discarded to avoid any startup transients in the data and typically a minimum of 10 repeated cycles were used 
for data analysis to ameliorate the effects of any measurement noise. A wide range of filters were tested to ensure 
the final results were not impacted by the filter cut-off frequency. Lateral-directional oscillatory tests were 
performed in a similar manner at six nominal values of angle-of-attack in the range 2° to 22°, two amplitudes of 5°, 
and 20°, and nine non-dimensional frequencies ranging from 0.05 to 0.35 or dimensional frequencies in the range of 
0.0098 Hz to 0.0685 Hz. 

Static data for normal force and pitching moment coefficients, over a ±10° range of sideslip, are shown in Figs. 2 
and 3. Normal force was not affected and fairly limited variation occurs in the pitching moment coefficient with 
changes in sideslip. Figures 4 and 5 show the same coefficients compared with wind tunnel measurements for the p 
= 0° alpha sweep. In comparison to wind tunnel results, the normal force compares fairly well. The slope compares 
well but the absolute levels are low. The distinct break in the wind tunnel results at stall are also not present in the 
water tunnel results. Due to the very low Reynolds number, the boundary layers are all laminar. The laminar 
boundary layers, coupled with the high aspect ratio design, mean that the wings will under perform, producing less 
lift than the higher Reynolds number wind tunnel data. The pitching moment results present levels that are offset 
from the wind tunnel data with the water tunnel data being more positive as the flow over the wings separates 
prematurely. 

Sample flow visualization results are shown for 0°, 10°, and 15° angle-of-attack at 0° sideslip. More cases and 
details can be obtained in the RHRC final report^^. Flow visualization obtained with dye ports is shown in Fig. 6 and 
with a movable wand in Figs. 7 and 8. These graphics show a well behaved flow field and normal stalled flow 
formations at 10° and 15° angle-of-attack where more turbulent flow appears near the rear wing root. 




Figure 1. Model of the NASA Hybrid Wing Body with sting built for static and 
dynamic testing in Rolling Hills Research Corporation water tunnel. 
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Figure 4. Comparison of water and wind 
tunnel normal force coefficient, 0={S°. 
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Figure 5. Comparison of water and wind 
tunnel pitching moment coefficient, 



Figure 6. Flow visualization at 0°, 10°, and 15° angle -of-attack and 0° sideslip. 


4 

American Institute of Aeronautics and Astronautics 





Figure 7. Flow visualization with dye wand over 
center body at 10° angle-of-attack and 0° sideslip. 


Figure 8. Flow visualization with dye wand off 
body at 10° angle-of-attack and 0° sideslip. 


IV. Model Structure with Unsteady Aerodynamics 

Model structures for the aerodynamic coefficients, , Cy , , and Q , represented by Q , where 

subscript a = N, m, Y, n, or I, were developed from a general form of the indicial model equations presented in Refs. 
5-7. A more general set of equations was developed for the wing-tail combinations allowing the wing and tail 
contributions to be treated separately in Refs 6 and 8. For the HWB configuration, however, a simpler model 
structure is used that allows the same form for both longitudinal and lateral -directional axes. In this case each 
lateral-directional coefficient can be written as 

Q (0 = Q (0) + 1 (t - T)P(r)dT + A J Q (? - T)p(r)dT + T)r(T)dT (la) 

0 0 0 

and the longitudinal coefficients can be written as 


Ca (0 = Ca (0) + j c {t - T)d(r)dT + j C (t - T)q{T)dr (lb) 

0 0 ^ 

In Eqs. (la) and (lb) coefficients , and are the indicial functions and 

Q(0) is the initial value of Q . From Ref. 7, two assumptions are adopted to simplify the model: (a) the effect of 
angular accelerations q , p , and r on any coefficient can be neglected and (b) the indicial functions in Eq. (1) can 
be expressed for the lateral-directional case as 

Cap (0 = (oo) - it) (2a) 


and for the longitudinal case as 


= (2b) 

where (t) and (t) are the deficiency functions and (qo) and (oo) are the rates of change with a 
and p , respectively, evaluated in steady flow conditions. 
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The simplified model, which takes into account changes with respect to steady state, has the lateral-directional 
form 

and the longitudinal form 

Cfl (0 = C (oo)«(0 + ^ c (oo)^(O - J F it - T)d{T)dT (3b) 

2v ^ Q 

To obtain a model appropriate for identification and with a limited number of parameters, the deficiency function is 
assumed to be a simple exponential function^ 

( 4 ) 

where ;rcan be cir or depending on axis being modeled. Models appropriate for an aircraft undergoing one degree- 
of-freedom forced oscillation in roll or yaw can be obtained using Eqs. (3) and (4). Considering one degree-of- 
freedom rolling motion in the tunnel 

= ( 5 ) 

where sideslip angle is related to the roll angle by the equation 

j3(t) - sin“^ (sin a{t) sin ^{t)) (6) 

Combining Eqs. (3-6), the aerodynamic models can be formulated as 

Q (0 = Cag («)A0 + Ca -aj p{T)dr (7) 

By introducing 

rj{t) = j P{r)dT (8) 

0 

and applying the Leibniz integral rule, the state space form of Eq. (7) can be written as 

f]{t) = -bdlit) + P{t) (9) 

Cait) = (<x.)A(0 + ^C„p {^)p{t)-aii(t) (10) 

Erom Eq. (7), a steady response can be obtained^ as 

Q (0 = sin((«0 + (j)A^ cos(cot) (11) 
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Where (j)A is the amplitude of roll oscillation, k is reduced frequency, and and are the in-phase and out-of- 

phase components, respectively. These components are related to the model parameters (aerodynamic derivatives) 
by the equations^ 


=Q^(oo)sin^Zo-a-^ 


rh^ 




sin^n 


( 12 ) 


Cap = Cap («>) - a ^ sin 


l + T^k^ 


( 13 ) 


For one degree-of-freedom yaw oscillations in the tunnel 

CAt) = CA¥it),r(t)] 

where yaw and sideslip angles are related by the equation 

j8(t) = sin“* (-cos a(t) sin 


( 14 ) 


( 15 ) 


The state space equations can be written as 




Ca it) = Ca^ i^)Pit) + — («)r(t) - ailit) 


( 16 ) 

( 17 ) 


and the corresponding in-phase and out-of-phase components' are expressed as 


Cap = Cap («) COS a^-a “o 


( 18 ) 


Ca, = Ca, (°c) + a cos «o 

1-i-rf^ 


( 19 ) 


For the longitudinal axis a similar derivation results in the following expressions. For one degree-of-freedom pitch 
oscillations in the tunnel 




( 20 ) 


where pitch and angle-of-attack are constrained by the physical test arrangement and thus related by the equality (for 
upright oscillations) 


a{t) = 0{t) 


( 21 ) 


The state space equations can be written as 


f]{t) = -b^r|{t) + d{t) 


( 22 ) 
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( 23 ) 


Ca (0 = Ca„ («)a(0 + — c {co)q(t) - arf{t) 

u 2\/ H 


and the corresponding in-phase and out-of-phase components^ are expressed as 


C =Q (oo)-a — 1- 


rh^ 


\ + Tik^ 


(24) 




l + T^k^ 


(25) 


V. Model Identification 

Two main steps in system identification are model structure determination and parameter estimation. Model 
structures presented in the last section have a form that allows parameters to be estimated using an output error 
method; however, some additional details on model order and degree of nonlinearity are required before parameter 
estimation is performed. Harmonic analysis is useful for that purpose by characterizing where, with respect to angle- 
of-attack and sideslip, unsteady behaviors occur and determining to what degree a nonlinear formulation is required. 
Harmonic analysis along with the coefficient of determination, can reveal these characteristics. is sensitive to 
noise so the analysis and interpretation of the coefficient can be made more difficult with low signal/noise data. 
Harmonic analysis is briefly described next; however, some discussion of the best practices for harmonic analysis 
and dealing with the effects of noise are presented more fully in Ref. 14. 


A. Harmonic Analysis 

A method of harmonic analysis^ ^ was applied to measured aerodynamic coefficients. A mathematical model for 
these coefficients is 


m m 

C^{t^ = Aq+ E Ay cos ( jcDt)-\- E Bj sin ( jcDt) a = l,n orY 


7=1 7=1 


(26) 


where Aq, Ay , and Bj are the Fourier coefficients. The analysis provides estimates of these coefficients, their 
standard errors, and the coefficient of determination, R^. For the model with linear aerodynamics and ^4 q = 0, the 
aerodynamic in-phase and out-of-phase components can be expressed in terms of the coefficients and . For the 
roll oscillation case the expressions are 


For the yaw oscillation case the expressions are 


c 


(27) 


C = 


k(pA 


(28) 


C =A 


(29) 






( 30 ) 
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where (5 is related to ^ and ^by Eqs. (6) and (15). For pitch oscillation the expressions are 


c =A 

«4 


( 31 ) 


C =A. 


ka. 


(32) 


The coefficient of determination, k, indicates the fraction of the variation in the measured data explained by the 
model and is defined as 


R^=\-SSe!SS, Q<R^<\ 


(33) 


where 




^^^ = Z[C_(/)-Q(0I 


(34) 


is the residual sum of squares and 


N 


;m2 


SS, = TVC,^ii)-C,{i)] 


(35) 


is the total sum of squares. (/),C^(/), and C^{i) are the measured, estimated, and mean values, respectively. 


B. Parameter Estimation 

Careful input design with a wide range of amplitudes and frequencies will provide sufficient information content 
for complete model identification. Parameter estimation for this study was accomplished using a variation on an 
output-error method^’^^’^^ in the time domain. In this approach the model equations are given by the state-space 
equations, noted in section IV, corresponding to the axis of interest. For the examples in this study, the longitudinal 
axis Eqs. (22) and (23) were applied. Measured time histories of the inputs, aor P , and outputs, Ca , for each of the 
frequencies tested are combined or “stacked” to ensure all the frequency content is included in the time domain 
estimation process. Equation integrations must be performed separately for each oscillation test but a single cost 
function is used to capture the overall model effectiveness. Application of the general output-error method for 
aircraft is explained and relevant software is provided in Ref. 19. 

VI. Modeling Results and Discussion 

Harmonic analysis was performed to test the aerodynamic model frequency response and to test the degree of 
nonlinearity for the HWB configuration. Sample results of harmonic analysis are shown in Fig. 9 for pitch axis 
oscillations. The values were obtained using the first harmonic and Eqs (31) and (32). The figure presents normal 
force and pitching moment out-of-phase and in-phase components as a function of angle-of-attack and frequency. 
The results indicate very little frequency spread (unsteady behavior) over the entire range of angle-of-attack for both 
normal-force and pitching -moment out-of-phase components and a limited spread for the normal-force in-phase 
component. The scales have been adjusted to magnify the very limited variation in magnitude at the different 
frequencies. 

Figure 10 shows the out-of-phase damping term for the HWB as a function of frequency, at four angles-of- 
attack, and a representative amplitude, = 10°. The lack of frequency dependence indicates that no unsteady 
behaviors occur for the conditions tested. In this case, without any frequency dependence present, investigators 
would normally test at lower frequencies to check for the characteristic “knee” in the curve of damping vs 
frequency. However, the extremely small values of at high frequency implies that this particular aircraft has 
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characteristically little damping in pitch and likely very limited unsteady responses even for lower frequencies. 
Future tests will verify this assumption regarding responses at lower frequencies. For a relative comparison, a 
conventional transport tested without the tail (only body-wing configuration) would present damping coefficients 
(out-of-phase coefficients) one to two orders of magnitude greater than presented in Fig. 9b, i.e., magnitudes in the 
range of 2 to 20. The scales in Fig. 10 are adjusted to a more representative magnitude of wing damping to 
demonstrate the frequency independence. 

Values of are consistently above 0.90 for the normal force coefficients indicating that a conventional linear 
aerodynamic model should be adequate for that coefficient. for the pitching moment coefficients are very low 
inidicating that a higher order model (nonlinear model) is required. Higher-order harmonics were found to improve 
the model fit. A 2”^-order harmonic model improved model fit such that R^ implied 70% of the variation was 

explained by the model and a 3^^ order model improved model fit such that R^ implied 81% of the variation was 
explained. Figure 1 1 shows the progressive improvement in R^ for the three different harmonic orders considered. A 
cubic nonlinearity is expected since the forced oscillation about cir^l0° with amplitude < 2 ^= 10 ° covers a range of the 
underlying static curve that is cubic in nature. 

The cubic characteristic found from harmonic analysis establishes a likely model structure for use during 
parameter estimation. Since the unsteady terms have been eliminated, the only source of nonlinearity must come 
from either static or steady-flow rotary dynamic terms. One option is to use stepwise regression, as used in Ref. 8, to 
select the best model structure. Another option, for low-order polynomials with only a few terms, is a trial-and-error 
approach that may be as fast and effective. The static model is available from measurements; so, estimating the 
static terms can be avoided by using the measured model directly. In that case, the identification problem can be 
simplified to estimating the best polynomial order for the damping terms. A progressive test of higher order terms 
can quickly reveal the model order required for the steady-flow damping. This approach resulted in a 3^^ order 
polynomial expression for the damping term as a function of angle-of-attack. Figure 12 shows measurements of the 
static pitching moment curve (blue) and forced-oscillation pitching -moment measurements (green) at each of the 8 
test frequencies. Each of the responses, shown as SCm versus Sa, have the nominal estimated static values subtracted 
out and represent displacements from the nominal test condition at ao=10° with oscillation amplitude ^^=10°. 
Estimated nominal static values from dynamic data are slightly different than the measured static values, hence the 
measured static curve does not pass exactly through 0 for Sa = 0. In addition the model-predicted responses (red) are 
shown where the steady-flow damping term has been estimated using a cubic polynomial representation. The 
steady-flow, nonlinear damping term in pitch was estimated as 

C„_^(oo;a) = -0.7155-3.7165a + 10.5296a^+111.1288c!:^ (36) 

A nonlinear polynomial for the static curve, within the range of oscillation, was also estimated for this case as 

(oo; cc) = -0.0060 + 0.01 36a + 0.3956a^ - 0.6384a^ (37) 

Overall, these models represent a very good fit to the measured responses at each frequency. Responses show the 
expected behavior of the dynamic response moving closer to the static curve as frequency is reduced. The estimation 
process is done including all frequencies simultaneously in order to ensure the best overall representation for time 
history responses. This was accomplished by “stacking” responses of each frequency in a modified output-error 
method discussed previously in section V.B. In this case, the R^ achieved by these models was 0.904 or in other 
words 90% of the variation was explained by the model over all the frequencies. Explaining 90% of the variation is 
an excellent result in any dynamic test environment. Poor measurement quality (signal/noise ratio) or degree of 
nonlinearity can limit the achievable R^ values. Polynomial coefficients and their standard errors are given in the 
table below. Mean values at the nominal or mean oscillation test point are subtracted from the measurements during 
output-error estimation so the first term, Aq , is not estimated for the static term C^(oo;ci') polynomial. The best 

estimate of the static polynomial term, Aq , is the static value of the pitching moment at ao. 
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Table: Polynomial coefficients and their standard errors for (oo;cif) and [oo;a) third-order models. 


Coefficients 

A) (A 

A (A 

A (A 

A (A 

(«>;«) 

N/A 

0.0136 (0.0007) 

0.3956 (0.0024) 

-0.638 (0.028) 

Cmg («;«) 

-0.7155 (0.0022) 

-3.717 (0.041) 

10.53 (0.21) 

111.13 (2.38) 




a, deg 

Figure 9a. 1^^-order harmonic analysis for 
Normal force coefficient, aA= 10°. 


a, deg 

Figure 9b. 1^^-order harmonic analysis for 
pitching moment coefficient, 10°. 
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Figure 12. Measured and computed pitching -moment responses at each test frequency, aA= 10°. 
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VII. Concluding Remarks 

Selected experimental results of a water tunnel study, at Rolling Hills Research Corporation, of the 0.7% scale 
NASA HWB-N2A configuration were presented in this paper. Static and dynamic tests about all three axes and flow 
visualization tests were performed using the HWB model. This work provides new dynamic data to support 
investigations at NASA Langley in developing advanced transport concepts and to support development of transport 
aircraft simulations in nonlinear or adverse flight regimes. 

Dynamic testing was performed using both single frequency and wide -band motions. Measured data were 
obtained over an extensive range of angles-of-attack, amplitudes, and frequencies. In this paper analysis was 
presented for single frequency forced oscillations in the longitudinal axis. The analysis found that HWB does not 
present significant unsteady behaviors but does show nonlinear characteristics for pitch damping and the static 
coefficients in the longitudinal axis. Cubic polynomials were required to model the steady flow damping term under 
large amplitude (10°) forced oscillations in pitch. 
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